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ABSTRACT 

Context. After many attempts over more than a decade, high-resolution spectroscopy has recently delivered its first detections of 
molecular absorption in exoplanet atmospheres, both in transmission and thermal emission spectra. Targeting the combined signal 
from individual lines in molecular bands, these measurements use variations in the planet radial velocity to disentangle the planet 
signal from telluric and stellar contaminants. 

Aims. In this paper we apply high resolution spectroscopy to probe molecular absorption in the day-side spectrum of the bright 
transiting hot Jupiter HD 189733b. 

Methods. We observed HD 189733b with the CRIRES high-resolution near-infrared spectograph on the Very Large Telescope during 
three nights, targeting possible absorption from carbon monoxide, water vapour, methane and carbon dioxide, at 2.0 and 2.3 fim. 
Results. We detect a 5-cr absorption signal from CO at a contrast level of ~4.5xl0~ 4 with respect to the stellar continuum, revealing 
the planet orbital radial velocity at 154+ 4 km s . This allows us to solve for the planet and stellar mass in a similar way as for stellar 
eclipsing binaries, resulting in M, = 0.846^'** Mq and M p = 1. 162+°°|| Afj up . No significant absorption is detected from H 2 0, C0 2 
or CH 4 and we determined upper limits on their line contrasts here. 

Conclusions. The detection of CO in the day-side spectrum of HD 189733b can be made consistent with the haze layer proposed to 
explain the optical to near-infrared transmission spectrum if the layer is optically thin at the normal incidence angles probed by our 
observations, or if the CO abundance is high enough for the CO absorption to originate from above the haze. Our non-detection of 
CO2 at 2.0 jim is not inconsistent with the deep CO2 absorption from low resolution NICMOS secondary eclipse data in the same 
wavelength range. If genuine, the absorption would be so strong that it blanks out any planet light completely in this wavelength 
range, leaving no high-resolution signal to be measured. 
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The extrasolar planet (exoplanet) transiting the K-star HD 
189733 is one of the most studied exoplanets to date, since 
the brightness of the system (V=7.75 mag) allows the de- 
tailed study of its atmosphere through transit and secondary 
eclipse measurements. Optical transmission measurements show 
a steady decrease of the effective planet radius from 0.3 to 
1 /im, which is thought to be due to Rayleigh scattering by 
a high altitude haze layer dPont et all 120081: ISing et al.1 [2009; 
iPont et al]|2012l) . The Rayleigh slope seems only interrupted by 
an absorption signal from sodium at 0.59 |/m (Redfiel d et alj 
2008; Huitson et al.l 1201 21) . There is a strong ungoing debate 
on whether this Rayleigh scattering signal continues in the 
near-infrared. Analysis of HST NICMOS spectroscopic ob- 
servations by (Sw ain et al.l 12008) show absorption signatures 
from water and methan e. However, reanalysis of the same 
data (Gibson e t~ai]|201 lb . narrow band NICMOS observations 



JSing et all [20091) . and HST data from WFC3 dGibson et all 
show no clear evidence for molecular features and suggest 



(Sing 

2012) 



* Based on observations collected at the European Southern 
Observatory (186.C-0289) 



that Rayleigh scattering dominates the transmission spectrum up 
to 2.5 yum. 

Numerous secondary eclips e me asurements (e.g. 



Deming et al. 2006; Grillm air et all 120081: ICharbonneau et al.l 
2008; ISwainet ail 120091) have placed constraints on the 
day-side temperature structur e and trace gas abundances (e.g . 
iMadhusudhan & Seageil2009tlLine et alJ2012tlLee et alJ2012b . 
In addition, phase curve observations from the Spitzer Space 
Telescope suggest that the hottest point in the planet atmosphere 
is shifted with respect to the sub-stellar poi nt, indicating strong 
longi tudinal circulation in the atmosphere (Knu tson et al.l l2009. 
|2012|) . However, a lso several s e conda ry eclipse measurements 
are controversial. Swa in et all (|2010) claim the presence of 
an enormous emission feature at 3.25 fim from ground-based 
observations, int erpreted as non-LTE e mission from methane. 
The same team (Wald mann et al.l [20121) confirmed this result 
with new data, in addition to sho wing simila r ly s trong emission 
features at 2.3 /mi. However, M andell et alj J2011I) see no sign 
of the L-band emission in NIRSPEC data, with upper limits 
40 times smaller than the claimed non-LTE methane signal, 
and suggest that imperfect correction for telluric water is the 
source of the feature claimed by Swain et al. In any case, the 
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unconventional observation and data analysis preformed by 
Swain et al. and Waldmann et al. needs to be validated using 
other techniques. 

The current measurements of the day-side of HD 189733b 
are taken at low spectral resolution at best, with most observa- 
tions being from broad-band photometry. The limited spectral 
resolution of the available measurements result in a large ambi- 
guity in atmospheric co mposition (e.g. IMadhusudhan & Seagerl 
2009; Line et al. 2012; Lee et al.ll2012l) and temperature struc- 
ture (Pont et afl l2012l) . 

Recently, high resolution spectroscopy has delivered its 
first detections of molecular absorption: carbon monoxide 
was detected in th e transmission spectrum of HD209458b 
dSnellen et al. 20101) and in the thermal daysid e spectrum of r 
Bootis b (jBrogi et al] |2012t iRodler et al1l2012l) . These observa- 
tions, conducted with the CRy ogenic high-resolut ion InfraRed 
Echelle Spectrograph (CRIRES iKaeufl et alj|200ll) on the Very 
Large Telescope (VLT) at a resolution of RxT00,000 target the 
combined signal from individual lines in the CO 2.3 //m molec- 
ular band. The large variation in the planet velocity is used to 
disentangle its spectral features from the telluric and stellar con- 
tamination. This also allows the determination of the orbital ve- 
locity of the planet, and therefore the orbital inclination and true 
planet mass in the case of non-transiting planets. R ecent results 
from thermal dayside spectroscopy of 51 Pegasi b (Brog i et alj 
2013) are somewhat ambiguous. Although the first two nights 
of data show a combined ~6-cr detection from carbon monoxide 
plus water, no signal is detected in the third night. Although the 
signal on the third night was expected to be slightly weaker due 
to the reduced quality of the data and contamination by stellar 
lines, a complete absence was not expected based on the first 
two nights. This would mean that either some yet unknown sys- 
tematics is affecting the data, or that time variability (by e.g. 
global weather patterns), plays a large role in the atmosphere of 
this object. 

High-resolution spectroscopy has been attempted before for 
HD 189733 b. Using two nig hts of Keck/NIRSP EC observations 
at a resolution of ~25,000, B arnes et al.l (1201 Oh identify a possi- 
ble planet signal using a deconvolution method. However, they 
could not make a strong clai m due t o the l ow significance of the 
result. Very recently, R odler et al.l (1201 3 1) revisited this dataset 
and claim a ~3-cr detection of CO lines in HD 189733b. 

Here, we present high-resolution CRIRES observarions of 
HD 189733 at 2.0 and 2.3 //m targeting absorption not only from 
CO, but also from CO2, H2O and CH4 in the planet's day-side 
spectrum. Section 2 describes the observations and data analysis 
and Section 3 the cross-correlation procedure used to combine 
the signal from the individual lines. Section 4 describes the re- 
sults, which are discussed in Section 5. 

2. Observations and data analysis 

We observed HD 18973 3b as part of a larg e programme (ESO 
programme 186.C-0289, Sne llen et al.ll201 1!) to probe molecular 
signatures in hot Jupiter atmospheres with CRIRES. CRIRES 
consists of four 1024x5 12-pixel Aladdin II detectors, separated 
by small gaps of roughly 100 pixels. The M ulti- Application 
Curva ture Adaptive Optics system (MACAO, lArsenault et al.l 
2003) was employed to maximise the throughput of the slit. 
Here, we use data taken during three nights, conducted at two 
wavelength settings, centred at 2.0 and 2.3 //m. The data were 
taken over a planet orbital phase-range near but not during sec- 
ondary eclipse to maximise the illumination of the planet and the 
planet's change in radial velocity (see Table Q~]for more details). 



Table 1. Details of the CRIRES observations, showing the ob- 
serving date, the targeted walength range, the planet orbital 
phase, the number of spectra after background subtraction, and 
the net integration time. 



Date 



/l(pm) 



Phase 



N 



(s) 



13 Jul 2011 
29 Jul 2011 
21 Aug 2011 



2.2875-2.3452 
1.9626-2.0045 
1.9626-2.0045 



0.3837-0.4801 
0.5727-0.6661 
0.3506-0.4377 



110 13200 
69 13800 
76 12160 



All observations used a 0.2" slit, resulting in a spectral resolu- 
tion of R« 100,000, and the star was nodded along the slit in an 
ABBA pattern to subtract the background. The spectra around 
2.0 fim contain possible planet lines of H2O and CO2, whereas 
the 2.3 fim region can be used to probe CO, H2O and CH4. 

On 13 July data was taken with an integration time of 4x15 
seconds per nodding position. On 29 July the integration time 
was 5x20 seconds and on 21 August 4x20 seconds. All three 
observations spanned roughly 5 hours each. 

The raw data was flat-fielded, background-subtracted, and 
optimally extracted using the standard CRIRES pipeline (version 
1.1 1.0). The data from the three nights, for each of the four de- 
tectors, were analysed separately. In each case, a 1024 x7V matrix 
was constructed from the N spectra. Bad columns in the matrix 
(corresponding to certain pixels in the spectra) were identified 
by taking the standard deviation of each column and finding >5- 
cr outliers. These bad columns (of the order of 20 per detector) 
were masked in the rest of the analysis. Further individual bad 
pixels were replaced by spline-interpolated values. 

The next step in the data reduction was to align all spectra 
onto a common wavelength grid. For this purpose we calculated 
a template telluric transmission spectrum with water, carbon 
dioxode and methane lines from HITRAN 2008 (Rothman et al. 
2009) using line-by-line calculations with a single model layer 
and pressures, temperatures and gas concentrations relevant to 
Paranal. Since HD 189733 A has significant CO lines in its spec- 
trum, we also calculated CO absorption lin es at the st ellar tem- 
perature using HITEMP absorption data (Roth man et al.| [2010) 
and scaled these relative to the telluric lines by comparisons with 
the data. We then Doppler-shifted the stellar CO lines accord- 
ing to the Earth's velocity, the known stellar systematic veloc- 
ity and its planet-induced radial velocity variation (Triau d let al.l 
2009), and combined them with the telluric spectrum. This was 
then correlated with each measured spectrum to find the best- 
fitting second-order polynomial correction on our initial guess of 
the wavelengths of each pixel, where the polynomial represents 
a wavelength-dependent shift. Finally, the spectra were spline- 
interpolated onto a common wavelength grid. 

Before we searched for molecular features in the data using 
the cross-correlation technique, we first needed to remove tel- 
luric and stellar signatures. The telluric lines are static in posi- 
tion, but vary in strength due to changes in the geometric airmass 
and variations in the water vapour content of the atmosphere. 
The stellar lines are invariable in time to a high level, but slightly 
change in position during the night due to variations in the radial 
velocity component of the Paranal observatory with respect to 
the star (mainly caused by the Earth's rotation). 

In this work we analyse the data i n a somewhat different w ay 
than previously dSnellen et al.ll201 It iBrogi et al.ll2012i 1201 3l) to 
explore new ways of reducing this type of data. We also per- 
formed the analysis with the previous method and found com- 
patible results (see Section 4.1). Instead of removing the telluric 
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contamination by fitting the columns in the matrix first with the 
geometric airmass and subsequently with the residual signals 
present in the strongest tellu ric lines, here we use singular value 
decompositions (SVD, e.g. Kalman 1996). For this we used the 
algorithm available in IDL. The SVD decomposes a matrix A 
into the following matrices: 



A = UWV T , 



(1) 



where U is the matrix of left singular vectors, W a (by def- 
inition) diagonal matrix of singular values and V the matrix of 
right singular vectors. In our case the matrix A represents the ar- 
ray of spectra on a single detector, with the number of pixels and 
number of spectra as dimensions. The right singular vectors are 
the eigenvectors and the left singular vectors, scaled by the sin- 
gular values, indicate how each of the eigenvectors are weighted 
to form each observed spectrum. A convenient property of the 
SVD is that matrix A can be approximated in an optimal way us- 
ing only the largest few sing ular values, and setting the smaller 
singular values to zero (e.g. An drews & Pattersonlll975l) . Here, 
we do the reverse and set the largest singular values to zero to 
remove the signals that are static with time while affecting the 
moving planet signal as little as possible. Fig. Q] shows an ex- 
ample of a matrix of spectra (detector 2 on the night of 13 July 
201 1) with progressive subtraction of the singular values for de- 
tector 2 at 2.3 fj.m. Examples of the left and right singular vectors 
are shown in Figs.|2]and[3] 

The left singular vectors corresponding to the first two sin- 
gular values of the 2.3 /vm data show clear correlations with the 
signal level in the continuum and the airmass for all four detec- 
tors (see Figs. |2]and|4j). Hence, the subtraction of the first two 
singu lar values is in eff ect very similar to the first steps used 
by ISnellen et all d2010t) and iBrogi et all ( 120121 120131) . Also the 
third left singular vector is very similar for all four detectors. 
Its corresponding eigenvector (Fig. |3) indicates it is related to 
shifts and/or broadening of the lines in the spectrum and its be- 
haviour with time is somewhat similar to parameters relevant to 
the adaptive optics performance (see Fig. |4). The only parameter 
in the headers that was similar to the fourth left singular vector 
is one of the temperature sensors in the instrument. The physical 
correspondances of these third and fourth left singular vectors is 
far from conclusive, with correlation coefficients of around 0.5 
between these vectors and the quantities plotted in Fig. [4] but 
it does highlight the diagnostic value of the S VDC for tracking 
potential systematic effects in the measurements. 

As a final step, we divide each column of the array of spectra 
by the square of the standard deviation. This prevents the cross- 
correlation analysis to be dominated by very noisy pixels. 



3. Correlation analysis 

We cross-correlated each reduced spectrum with model spec- 
tra of exoplanet atmospheres with Doppler shifts ranging from 
-200 to +200 km/s relative to the expected planet ra dial velocity- 
Orbita l parameters of the system were taken from [Triaud et al J 
(2009). In order to scan the parameter space and check for po- 
tential spurious correlation signals, we also performed our cal- 
culations for a wide range of amplitudes of the radial velocity 
of the planet, K p (50-200 km s ). We assumed a circular orbit 
throughout the paper, which is reasonable given that the maxi- 
mum offset in plane t radial ve locity for its possible small eccen- 
tricity (see iTriaud et al.ll2009l) is smaller than the velocity reso- 
lution of CRIRES. 
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Fig. 1. Example of stellar and telluric signal removal: the array 
of spectra (in ADU per second) for detector 2 at 2.3 //m on 13 
July 201 1 after setting the first n singular values to zero, with n 
being indicated above the panels. 



The model spectra are the result of line-by-line calculations 
using a range of parameterised pressure-temper ature structures 
(see b e low), H2-H 2 collision-induced absorption (Bor vsow et al.l 
2001; iBorvsowl |2002) and absorption from a single trace 
gas. Line data for the trace gases are taken from HITEMP 
2010 for H 2 Q and CO dRothmanet alj 12010b. HITRAN 2008 
for CHa JRothman et al] |2009j) and HITEMP 1995 for C0 2 
(Rothman et aljfl995l) . A Voigt line-shape is used for calcula- 
tion of the absorption. The atmospheric temperatures at pres- 
sures higher than 0.1 bar, where we take a temperature of 1350 
K, are made to follow t he bes t fitting temperature profiles of 
Mad husudhan & Seagerl J2009I) and the temperature at a lower 
pressure p\ is varied from 500 to 1500 K in steps of 500 K. 
As such, we also included the option of weak thermal inver- 
sions at high altitudes. Between 0.1 and p\ bar a constant lapse 
rate (i.e. the rate of temperature change with log pressure) is as- 
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Fig. 2. The first four left singular vectors for all four detectors 
(black - detector 1; red - detector 2; green - detector 3; blue 
- detector 4) at 2.3 jum on 13 July 2011, divided by their me- 
dian value. The corresponding eigenvalues for detector 2 are also 
given to indicate their relative strengths. These vectors can be 
compared to the parameters in Fig. [4] 
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Fig. 4. Instrument and environmental parameters possibly re- 
lated to the left singular vectors plotted in Fig. [2] The median 
count is for detector 2 at 2.3 //m on 13 July 2011, AO Strehl 
indicates the Strehl ratio from the adaptive optics system in per- 
cent and S V temp is one of the many temperature sensors on the 
instrument, located at the slit viewer. 
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Fig. 3. The first four right singular vectors, or eigenvectors, for 
detector 2 at 2.3 /tm on 13 July 201 1. 



sumed and p\ is varied from 10-10 in multiplication steps 
of 10 . At pressures lower than p\ a uniform temperature is 
assumed. Volume mixing ratios of the gases ranged from 10 -6 - 
10~ 3 in multiplication steps of 10 05 . 

The correlation analysis was performed for each of the four 
detectors separately. To combine the correlation signals obtained 
from different spectra and detectors we first determined how to 
weigh their contributions. We did this by inserting an artificial 
planet signal in the wavelength-aligned data for the system ve- 
locity and K p under consideration, and perform the same data 
reduction as the original data. If the data with inserted signal 
gives a large increase in correlation value for a particular spec- 
trum, this spectrum is thus expected to give a high signal also 
in the original data. Hence, the increase of the correlation when 
inserting the artificial planet signal can be used as a weighing 
factor to favour spectra or detectors for which a larger sign al is 
expected. The correlation values for each spectrum and detector 



are multiplied by this increase of correlation, and added together 
to form the total correlation signal. 

In the correlation analysis for CO, there are still some resid- 
uals originating from the stellar spectrum. This is not surprising, 
since the star's radial velocity slightly changes with time due to 
the Earth's rotation and the star has significant CO lines in its 
spectrum. Fortunately, the velocities of the planet are at least 20 
km/s offset from that of the star for the phases observed, so the 
strongest stellar residuals do not affect the correlation signal of 
the planet. When combining the correlation signals of the differ- 
ent spectra, we mask the velocities that are within 5 km/s of the 
stellar velocity after the cross-correlation. Note that this does not 
significantly affect our cross-correlation signals at the expected 
planet velocity, but only within the region of the expected stellar 
signals. 

As more singular values are subtracted from the data, parts 
of the planet signal may also be removed. Hence, there is an op- 
timum number of subtractions that will lead to the best recovery 
of the planet signal. We tested this effect by again inserting arti- 
ficial planet signals and evaluating the correlation signal for the 
different numbers of reduced singular values. Fig. |5] illustrates 
how the planet signal in a single detector first becomes apparent 
when the stellar and telluric signals are removed, and then slowly 
disappears as more and more singular values are removed. The 
slow drop of correlation with more subtracted singular values 
rises from the fact that there is only very little power in each 
of these higher singular values. The correlation as a function of 
number of subtracted singular values is actually very similar to 
the correlation of the data without inserted signal, but the corre- 
lation values are higher. For each detector, we chose the number 
of singular values to subtract according to the highest correla- 
tion of the inserted signal with respect to the correlation noise, 
when combined over the entire night. Fig.|5]also illustrates what 
happens to the signal if negative planet signals are inserted. Such 
a negative signal would correspond to spectra with lines seen in 
emission instead of absorption. 
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Fig. 5. Total correlation at the planet position as a function of 
number of subtracted singular values for detector 2 at 2.3 fim 
on 13 July 2011 for the planet model giving the highest corre- 
lation. Different lines indicate different multiplication factors of 
the inserted planet signal (see text), as indicated by the labels. 



4. Results 

Our correlation analysis yields a detection of carbon monoxide 
at 2.3 jt/m, which is presented below. No significant signals are 
detcted for the other bands and molecules, of which the implica- 
tions are discussed in Section 4.2. 



4.1. Detection of CO 

The cross-correlation assuming different system velocities and 
values of the planet radial velocity K p for our best CO template 
are shown in Fig. [6] The best CO template itself is shown in 
Fig. [7] A positive correlation value (black) means in this case 
that the CO lines are apparent in absorption, as in Fig. [7] indi- 
cating a non-inverted temperature profile over the probed pres- 
sure range. The correlation shows a 5.0-cr maximum at the ex- 
pected system velocity, within the range of K p expected from the 
measured radial vel ocity of the star an d its estimated mass from 
spectral modelling dTriaud et al.ll2.009b . The significance was es- 
timated by the peak cross-correlation value divided by the noise 
in the cross-correlation values in Fig. 6 at system velocities away 
from the planet signal. The signal is the combination of all four 
detectors, although detector 1 does not show a clear signal. Note 
however that detector 1 has a low weight, since artificially in- 
jected planet spectra give only very faint retrieved signals. This 
is consistent with the CO spectrum, which shows less lines in 
the wavelength range covered by detector 1 . 

The observed CO signal is detected at K p - 154^1 km s _1 
and V sys = -2 km s" 1 , both consistent with the lite rature values 
for th e orbital and stellar parameters of the system (fTriaud et alj 
2009). We estimated the statistical significance of the signal 
in several ways. First of all, we tested how the distribution of 
cross-correlation values as calculated for the 110 spectra com- 
pare to a Gaussian distribution. A histogram of the correlation 
values as obtained from detectors 2-4 is shown in Fig. [8] The 
interval of the system velocities over shich the histrogram was 
plotted was chosen such that the area of the stellar signal was 
avoided and a large sample was obtainted. The best fit of the 
histogram to a Gaussian distribution is overplotted, showing a 
very good agreement. A further analysis of the quantiles of this 
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Fig. 6. Combined cross-correlation values for all four detectors 
as a function of system velocity and K p . The dashed line indi- 
cates where the stellar CO lines are expected to leave the largest 
contamination signal. The area in between the dashed lines can 
also be contaminated by stellar residuals. The dotted lines indi- 
cate the system velocity of HD 189733 and the peak position of 
K„. 
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Fig. 7. Template spectrum that gives the best correlation value. 
Note that we are not sensitive to the continuum level of the spec- 
trum. For the sake of clarity, stellar lines are ignored in the cal- 
culation of the planet-to-star contrast F p /F s . Brackets indicate 
the wavelength ranges covered by the four detectors. 



distribution (Fig. |9]l shows it exhibits no significant deviation 
from Gaussianity. Only detector 1 shows some non-Gaussian 
behaviour using the optimal number of singular values, which 
improved when more singular values were subtracted. Since we 
did not find a significant signal in this detector and it is weighted 
very lightly compared to the other detectors, and because the ex- 
tra singular values effectively only remove very little of the data, 
the additional subtraction of the singular values had no impact 
on the final result. 
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Fig. 8. Histogram of correlation values outside the expected 
planet velocity (+10-400 km s _1 with respect to the planet) for 
detectors 2-4 combined. Error bars are the suare root of the n um- 
ber of occurrences in each bin. The dotted line shows a Gaussian 
fit. 
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Fig. 9. Quatiles of a normal distribution versus the measured 
quantiles of Fig. [8] The gray line indicates a perfect normal dis- 
tribution. The two agree well up till the 4-cr level, which is de- 
termined by the number of available sampling points. 



An additional test we performed was the Welch T-test (Welch 

1947), comparing the planet signal distribution with the correla- 
tion distribution away from the expected planet velocities (see 
Fig. [Tol l . Unfortunately, we have far less spectra available com- 
pared to Brogi et al.l(l2012il2013l) . who performed the same test, 
so the sampling of the correlations at the planet velocity is worse 
here. Nevertheless, we can reject the null-hypothesis that the two 
samples are drawn from the same distribution at a 4.2-cr confi- 
dence level. 

Our measurements are not sensitive to the continuum level 
of the planet spectrum, since we remove any broadband signal in 
our data analysis. We are only sensitive to the depths of the nar- 
row absorption lines, which are the parts of the planet spectrum 
that show the Doppler shift. We determined the line contrast (the 
depth of the deepest lines with respect to the continuum, divided 
by the stellar flux) by inserting the model that gives the strongest 
signal times different (negative) scaling factors (see Section 3) 
until a significance value of 0+1 <x was reached at the planet 



Fig. 10. Histogram of correlation values away from the expected 
planet velocity (+10-400 km s -1 with respect to the planet), and 
those at the planet velocity (gray line, scaled), for detectors 2-4 
combined. 



velocity for all four detectors combined. The line contrast for 
the deepest CO lines determined in this way is (4.5±0.9)xl0~ 4 , 
which corresponds roughly to the signal from the template spec- 
trum (see Fig. 13. This is the value for the model spectra and does 
not include the convolution to the instrument resolution. The line 
contrast including convolution gives a value of (3.4±0.7)xl0~ 4 . 
This line contrast can be interpreted as the difference of the line 
cores with respect to the broadband continuum and can be used 
as an additional constraint on the temperature profile and CO 
concentration (see the discussion in Section 5). The continuum 
level will be close to the broadband secondary eclipse depth 
around 2.3 fim for this transiting planet. 

Using the data analysis method we pr esented previously 
(ISnellenetal.H2010t iBrod et al.ll20Tl 120131) we found a 4.8-cr 
detection of the CO signal, which is consistent with the results 
presented here. Hence, either method can be used equally well, 
although the method presented here is perhaps less ad hoc and 
allows for a more objective treatment of the data. 

4.2. Upper limits for C0 2 ,H 2 and CH A 

Besides CO, we also searched for signals from H2O and CH4 
in the 2.3 //m data, and CO2 and H2O in the 2.0 //m data, but 
did not find a significant signal. We determined upper limits of 
the line contrast by inserting inverted models until a significance 
value of -3<r was reached at the planet velocity for all four detec- 
tors combined. We used the same number of singular values as 
for the CO analysis. The resulting upper limits of the maximum 
model (unconvolved) line contrasts is 2.8xlO -3 for H2O and 
8.3xl0~ 4 for CH4. The latter is comparable to the value of the 
medium-resolution secondary eclipse depth of 8xl0~ 4 , i.e. the 
total planet-to-star flux contrast at th ese wavelengths, a s mea- 
sured with NICMOS on the HST by I Swain et all J2009I) . Since 
the continuum level would be close to the medium-resolution 
planet/star flux contrast, this would imply that the flux in the 
CH4 lines would have to drop to zero for them to be detectable at 
the determined upper limit of the line depth. The contrast of the 
best-fitting CO lines is ~ 4.5xl0 -4 , so CH4 would have to be sig- 
nificantly more abundant than CO for this to happen, given their 
similar line strengths in this region. Also, the upper atmosphere 
of HD 189733b would need to be extremely cold (of the order 
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of a few hundred Kelvin). Neither of them is very li kely (see 
e.g. iMadhusudhan & Seagerl 120091: iMoseset alj|2013l) . Hence, 
the upper limit on the line contrast of CH4 does not constrain 
the abundance in a meaningful way if the measured medium- 
resolution planet-to-star contrast is correct. 

Although we have collected two nights of data at 2.0 //m, 
we were less sensitive to atmospheric signals in this wavelength 
range due to a combination of more telluric absorption and a 
smaller contrast between the planet and the star. Fig. [TT] shows 
an example spectrum of the 2.0 /vm data. The spectra on detec- 
tors 1 and 4 are dominated by telluric CO2 lines, greatly reducing 
the potential for retrieving exoplanet signals. In detectors 2 and 3 
also water lines are evident. In general, the telluric transmission 
is significantly lower than around 2.3 /urn. For this wavelength 
setting we ignored detector 4 due to the much lower count level. 
The SVD analysis could reduce the systematic signals down to 
similar levels as in Fig.Q] Again, models with varying tempera- 
ture profiles and volume mixing ratios of H2O and CO2 were in- 
serted with K p - 154 kms~' tofind3-crupperfimitsof 8.7xl0~ 4 
forH 2 and 1.8X10 -3 forC0 2 . 
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Fig. 11. Example spectrum of the data at 2.0 jt/m. 



4.3. Planet system parameters 

Since we now have a measure of the velocity and the orbital in- 
clination of both the star and the planet around their common 
barycentre, model-independent estimates of their masses can be 
obtained in the same way as has been don e for stellar eclips - 
ing binaries for more than a century ( e.g. iTorres et al.l 12 010). 
Using the values and uncertainties from lTriaud et al.l(l2009[ see 
Table 2) for the amplitude of the stellar radial velocity vari- 
ations, planet orbital inclination and period, we determine a 



stellar mass of M, 



0.846 ^j;;^ M© and a planet mass of 
M p — 1.162^1"^ Mj up . The stellar mass is consistent with that 
found bv lTriaudetaH d2009l) (M s = 0.823+°;°g M Q and M p = 
I.I38+0025 ^Jup)- In principle, if these CRIRES observations are 
repeated on both sides of the secondary eclipse (phases higher 
and lower than 0.5) and also at somewhat higher planet radial 
velocities (phases closer to 0.25 or 0.75), the uncertainty in stel- 
lar mass can be reduced to a level of ~1%, This issue is further 
discussed in Section 5. 



Table 2. Stellar and orbital parameters taken from iTriaud et al.l 
(2009) and derived from our measurements. 



Parameter (unit) 



Value 



Assumed parameters (Triaud et al., 2009): 

^(ms -1 ) 

P (days) 

in 

Measured parameters: 
K p (km s- 1 ) 
Derived parameters: 
M s (Mo) 

M P (M luv ) 



201 96 +1 ° 7 
z.vi.7u _ 063 

? ? 1 8S7^ 1 9+0.00000036 
Z..Z.10J / Jli_„ 00000076 



85.508 



-0.05 



154+| 



0.846*"» 

1 162 + ™*s 

1,10 -0.039 



5. Discussion and conclusions 

Using high resolution spectra from CRIRES on the VET at 
orbital phases between 0.38-0.48, we detect carbon monox- 
ide in the day-side atmosphere of HD 189733b at a 5-cr 
significance. Fro m Spitzer data the presence of C O was al- 
ready suspected (lBarmanll20"08l: iDesert et all 120091) , although 
more wider scans of the parameter space revealed no use- 
ful constrai nts on CO based on the Sp itzer secondary eclipse 
data alone JMadhusudhan & Seage r 2009). HST/NICMOS sec- 
ondary eclipse dat a do show large CO volume mixing ra- 



iios (~10- 4 -lQ- 2 ) dSwain et al.l 120091: IMadhusudhan & S eager 



l2009t iLee et al.l 120121) , although past NICM OS results hav e 
been received w i th som e scepticism (e.g. Gib son et al.ll201 ll) . 
IWaldmann et al.l J2012I) claim to detect a strong planet emis- 
sion feature at 2.3 fim using low-resolution groundbased spec- 
troscopy, similar to the controversial signal found at 3.2 fim by 
ISwain et al] d2010l) . and also suggest it could be non-LTE emis- 
sion. Due to the fact that no proper physical model exists yet for 
this emission, we could not perform a detailed search for this 
emission in our data. However, no obvious high-resolution sig- 
nal, expected from the ~2xl0~ 3 low resolution signal seen by 
Waldmann et al., is present in our data. Figure 1 indicates that 
single lines cannot be seen above the 1 % level in our data. 

The CRIRES measurements are not sensitive to the contin- 
uum level of the planet spectrum, since we remove any broad- 
band signal in our data analysis. Instead, we determine an uncon- 
volved line contrast (the depth of the deepest lines with respect 
to the continuum, divided by the stellar flux) of (4.5+0.9)xl0~ 4 
for the strongest CO lines. This line contrast indicates the differ- 
ence between the thermal flux coming from the continuum level 
at high pressure and the flux from the line cores, which probe 
very low pressures (see Fig. [T2V The line contrast thus depends 
on the difference in temperatures at the pressures of the contin- 
uum and the line cores. Already it is clear that there is a large 
degeneracy in the temperature structure of the atmosphere that 
can reproduce a certain line contrast, since different sets of tem- 
peratures can have the same difference in flux. Furthermore, the 
abundance of CO adds another ambiguity, since increasing the 
CO abundance will decrease the pressure that is probed by the 
line cores. It is therefore uncertain at what pressure the flux of 
the line cores originates. This also means we cannot determine 
CO abundances from line contrasts alone. Fig.[T3lillustrates this 
point. It shows the maximum line contrasts that were obtained 
in our model grid (see Section 3). Note that in this grid we kept 
temperatures at pressures above 0. 1 bar fixed, meaning we ba- 
sically fixed the continuum level, already reducing part of the 
possible degeneracies. Furthermore, our grid only covered a se- 
lected portion of the parameter space. A lower limit on the line 
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contrast comes from an isothermal atmosphere, which would not 
show any lines (zero contrast). Fig.[T3lshows that there is prac- 
tically no dependence of the range of possible line contrasts on 
the CO volume mixing ratio. Hence, more information is needed 
before the CO abundance can be determined from these kinds 
of high spectral resolution measurements, most notably the tem- 
perature profiles covering pressures from the continuum level 
to the low pressures of the line cores. Constraints on the tem- 
perature profile can be obtained by measuring the thermal flux 
(or secondary eclipse depths) at a range of wavelengths, prob- 
ing a range of pressures due to differences in gas absorption 



Madhusudhan & Seager 2009; Li ne et all 120121: iLee et all 



|2012|) . However, lower resolution measurements do not probe 
the low pressures of the line cores well. Also, there is still a 
degeneracy with absolute gas abundances, since this determines 
the pressure levels that are probed. Ideally, one would have a 
single gas for which we know (or can estimate reliably) the ab- 
solute abundance, and which probes a wide range of pressures 
from absorption lines or bands of different strengths. From mea- 
surements of this gas at different wavelengths the temperature 
structure can then be determined, if the variation of the gas 
abundance with pressure is known. CO2 plays this role in the 
Venus atmosphere and CH4 in the atmospheres of Jupiter, Saturn 
and Titan. Obtaining absolute gas abundances for giant exo- 
planet atmospheres is possible using transit (transmission) mea- 
surements if the H2 Ray leigh scattering slope is measured (e.g. 
Benn eke & Seagerf 2012). although this signature might be diffi- 
cult to distinguish from that of small haze particles. In summary, 
obtaining absolute gas abundances from line contrasts such as 
the ones we present here is difficult, but the line contrasts can be 
used as an additional constraint to reduce the uncertainty in re- 
trievals. However, if multiple gases would be detected in a single 
high resolution observation, this would make it possible to de- 
termine the relative abundances of these gases, since they would 
both be reliant on the same temperature profile with the same 
continuum level at the same time. There might still be some de- 
generacy caused by the line cores of the two gases probing two 
different ranges of pressures, but these pressure ranges are likely 
to overlap. 
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Fig. 12. Normalised contribution functions for a model spectrum 
with CO at a volume mixing ratio of 10~ 4 . These show where in 
the atmosphere the thermal flux is emitted for each wavelength. 



In comparing our results with the claimed detection by 
Rodl er et al.l ([2013), we note that the planet radial velocity mea- 
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Fig. 13. The maximum line contrasts (solid line) as a function of 
CO volume mixing ratio in our model grid (see Section 3). The 
minimum line contrast in a planetary atmosphere is zero for an 
isothermal atmosphere. Dashed lines show our determined limits 
on the CO line contrast. 



surements ar e con sistent within their reported l-<x error bars. 
iRodler et ail (120131) find a CO line contrast of 1 .8 x 1 0~ 4 at a spec- 
tral resolution of R=25,000 with NIRSPEC at Keck. Assuming 
that the lines are unresolved, this translates to a line contrast of 
-7x10 4 at the CRIRES resolution, which is more than 50% 
higher than what we find in our study. This may suggest that 
the reliability of the Keck signal is somewhat o verestimated. 
Howev er, the error bars on the line contrast of IRodler et al.l 
(2013) will probably be large, although these are not given in 
their paper. R odler et aLl (2013) suggest a possible low CO abun- 
dance based on their analysis of the Keck data. Our line con- 
trast is potentially even lower, but unfortunately we cannot yet 
draw conclusions on the CO abundance based on the line con- 
trast alone (see the discussion above). Taking their 3-cr detection 
using 13 hours of NIRSPEC observations at face value, the 5- 
cr detection presented here implies that CRIRES is a factor 7 
more efficient. We attribute this to the high spectral resolution 
and stability of CRIRES, and the well-performing adaptive op- 
tics system on the VET. 

Based on the available secondary eclipse data, retrieval 
methods favour temperature profiles that decrease in tempera- 
ture as pressure decreases, i.e. a te mperature profile without an 
inversion layer (Madhusudhan & Seager 2009; Line et al. 2012; 
ILee et al]|2012l) . However, in recent work Pont et al.l (2012) ar- 
gue that, based on their transmission measurements and im- 
proved analysis of secondary eclipse data, HD 189733b could 
also have temperatures that are increasing with decreasing pres- 
sure. This would then be the result of an optically thick haze 
layer, for which they see the signature in the transit measure- 
ments. However, our detection of CO shows that the tempera- 
ture profile is not inverted at the pressures probed by the CO 
lines (10 5 -10 3 bar for a CO volume mixing ratio of 10~ 4 , as 
seen in Fig. Q~2] and roughly an order of magnitude lower pres- 
sures for a CO volume mixing ratio of 10~ 3 ) and that the CO 
lines are not obscured by an optically thick haze layer. In prin- 
ciple it is possible that we are only seeing the cores of the CO 
lines sticking out above the haze in a higher part of the atmo- 
sphere where the temperatures decrease again with altitude, but 
this would require large C O volume mixing ratios (> 10~ 4 ) if 
the temperature profile of iPont et all (1201 2|) . which shows an 
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inversion down to at least 0.1 mbar, is correct. Such a high 
CO abundance with absorption features above the haze layer 
would then also be visible around 5 fim, where CO absorp- 
tion is stronger than at 2.3 /urn. Such a CO feature might in- 
deed be visible in Spitzer data dKnutson et al.l 120121) . An op- 
tically thick haze layer would also require a cold atmospheric 
la yer above the haze layer, which is contrary to the conclusions 
of Huitson et al. (2 0120. Hence, we conclude that the haze layer 
seen bv lPontet all d2008l) : ISing et all d2009l) : iPont et all d2012l) 
is either optically thin when seen at normal incidence angles and 
that consequently there is no temperature inversion in the mid- 
dle atmosphere, or that the haze layer is optically thick, but low 
enough for CO lines to be visible. 

We could derive only upper limits for the line contrasts of 
H 2 0, C0 2 and CH 4 : 8.7xl0~ 4 for H 2 and 1.8xl0" 3 for C0 2 
at 2.0 fim and 2.8xl0~ 3 for H 2 and 8.3xl0~ 4 for CH 4 at 2.3 
fim. Future observations, also at more optimal wavelengths, will 
be needed to more effectively search for these gases and assess 
their relative abundances. Less telluric absorption and a higher 
planet-to-star contrast would be benificial in these observations. 
However, it is not straightforward to identify an optimal wave- 
length for these observations. For instance, the planet-to-star 
contrast is better at longer wavelengths, but there the sky back- 
ground increases. Also, there is a trade-off between targeting 
strong lines in the planet spectrum, while avoiding strong lines in 
the telluric spectrum. A sensitivity study covering the CRIRES 
wavelength range, using realistic instrument performance and 
telluric absorption spectra, is needed to determine the optimal 
wavelengths for detection of these gases using high resolution 
spectroscopy. Note that the optimal observing wavelengths may 
differ for different targets. 



154+ 4 km s- 1 
we derive a stellar mass of M s = 0.846 + nn*n Mq and a planet 



From the measured radial velocity of K p 



-0.049 

mass of M p - 1 J62*!M5g Mj up , independent of stellar spectral 
modelling. These values are consistent with previously known 
values and their errors are somewhat larger. However, combin- 
ing these results with possible future observations at phases 
larger than 0.5, as well as phases closer to 0.25 and 0.75, would 
help significantly in reducing the erro r on K D , and hen ce in 
the masses of the star and planet (e.g. iBrogiet al]|2013l) . An 
error in K p of 1 km s would result in an error on the stel- 
lar mass of 0.017 Mq, which is more precise than the cur- 
rent value derived from spectral modelling. The ability to treat 
transiting planet systems in the same manner as double-lined 
eclipsing binaries (DLEBs) allows a model-independent mea- 
surement of the stellar and planet masses and radii, which is 
particularly useful for planets with M-dwarf host stars. Precise 
dynamical measurements of the masses and radii of low-mass 
stars in DLEBs are not well-reproduced by stellar evolutions 
models and have radii that are up to 15% la rger than predicted 
dLopez-Morales & Ribasl 120051: iRibasI 120061) . Strong magnetic 
fields in the short-period low-mass DLEBs are thought to be re- 
sponsible for the rad ius inflation (see e.g. IChabrier et al.l l2007; 
Mor ales et al.l 12010). Even in higher mass G- and K-dwarfs, 
strong magneti c activity in the star s results in disagreement with 
stellar models (M orales et al.ll2009l) . However, in wider separa- 
tion M-dwarf DLEBs with magnetically inactive low-mass stars , 
the radii are still inflated (Ilrwinetal.ll201li iDovleet al.ll201ll) . 
with no real convergence towards stellar models at longer pe- 
riods (Bir kbv et al.ll2012l) . Obtaining K p for planets transiting 
M-dwarf s removes the uncertainty in stellar models and gives a 
direct and accurate characterisation of the star-planet system. 
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